Subwavelength focusing inside an open disordered medium 
by time reversal at a single point antenna 
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We study theoretically light focusing at subwavelength scale inside a disordered strongly scattering 
open medium. We show that broadband time reversal at a single point antenna, in conjunction with 
near-field interactions and multiple scattering, produces spatial focusing with a quality comparable 
to that obtained in an ideal closed cavity. This gives new perspectives for super-resolved optical 
imaging and coherent control of single nanosources or absorbers in complex media. 

PACS numbers: 42.25.Dd 



The ability to focus light in a small region of space 
underlies many optical imaging techniques, including the 
most recent improvements in optical microscopy Fo- 
cusing beyond the diffraction limit is a major objective, 
for the improvement of spatial resolution. One possi- 
ble approach is offered by scanning near-field optical mi- 
croscopy, in which a sharp tip produces a subwavelength 
spot that can be scanned in the near field of the medium 
under investigation Q. Recently, the combination of 
wavefront shaping and scattering by a structured envi- 
ronment has been recognized as an alternative way to 
produce subwavelength focal spots in optics An 
efficient control of wave propagation in complex media, 
including focusing, can be achieved by time reversal. The 
technique is well established in acoustics ML an d has been 
extended to seismology [§], microwaves |i| , and to op- 
tics [l(J 11 1- The field of subwavelength light focusing 
encompasses a number of fields beyond imaging, such as 
coherent control of single emitters or absorbers in com- 
plex media [l2| or addressing in integrated optics a ni- 
ii. 

A typical time-reversal focusing experiment consists in 
two steps. In the direct situation, the field radiated by a 
(point) source in an arbitrary medium is recorded on an 
array of detectors (often called a Time-Reversal Mirror 
or TRM). In the time-reversed situation, the detectors 
become sources sending back the time-reversed sequence 
of the recorded field. The result is a time-reversed field 
that focuses towards the source location. Time reversal 
below the diffraction limit has been demonstrated for the 
first time in acoustics, using an acoustic sink (i.e. an ac- 
tive time- reversed source) placed at the focal point [lij ]. 
In a structured environment, subwavelength focusing can 
be achieved even without creating a sink due to the pres- 



ence of scatterers near the focal spot [17[. Since recent 
technical progress has made possible the use of time re- 
versal concepts for spatial (and temporal) focusing of 
light in complex media [l8|, it seems that some funda- 
mental questions need to be examined: In particular, 
the conditions to obtain subwavelength focusing in an 
open disordered medium, and the possibility to perform 
a time-reversal experiment with a single antenna by tak- 



ing advantage of multiple scattering. 

In this Letter, we address these questions theoretically 
based on numerical experiments. We show that subwave- 
length focusing in a strongly scattering disordered open 
system is feasible by time reversal at a single point an- 
tenna. We put forward the crucial role of the spectral 
bandwidth, in conjunction with near-field interactions 
and multiple scattering, to reach focusing performances 
comparable to that produced in an ideal closed cavity, 
and that could not be achieved by monochromatic phase 
conjugation (l9| . 

We consider a two-dimensional cluster of cylindri- 
cal scatterers randomly distributed inside a cylindri- 
cal region of radius R = 1.91 pm. One specific con- 
figuration of the cluster is used throughout the paper 
[Fig. [T] (a)]. A minimum distance d m i n = 10 nm is 
forced between scatterers to avoid overlapping. The 
scatterers are described by their electric polarizability 
a (uj) = -2Tc\l [cj 2 (w - ujq + iT/2)] where uj = 3.14 x 
10 15 rad.s _ is the resonance frequency, T = 10 14 s _1 is 
the linewidth and Co is the speed of light in vacuum. This 
corresponds to the Transverse Electric (TE) polarizabil- 
ity of a resonant non-absorbing 2D-point scatterer, with 
a quality factor Q = ujq/T = 31.4. 

The time reversal process is assumed broadband, in 
a spectral interval [f2 — A, O + A], with CI the central 
frequency and 2A the bandwidth. For the numeri- 
cal simulations, we have chosen Cl = ljq — 5T and 
A = 10T. This corresponds to a reference wavelength 
A = 714 nm. The number of scatterers is TV = 11205 
such that the average distance between scatterers is (d) = 
32 nm. At the reference frequency CI, the Boltzmann 
mean-free path is £b {CI) = [pa s (CI)] 1 = 227 nm with 
a s (oj) = (CI / c ) 3 / 4 \a (Cl)\ 2 the scattering cross-section 
and p the density of scatterers. The optical thickness is 
b = 2R/£b (CI) = 16.8, large enough for multiple scatter- 
ing to take place and (d) /A = 0.045, small enough for 
near-field interactions to occur. The scattering strength 
is (Cl/c )e B (CI) = 2. 

In the direct situation, the 2D system is illuminated 
using a point source [dipole moment po (cu) = cte] po- 
larized along the scatterers (TE modes), such that the 
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FIG. 1: System composed by a 2D cluster of cylindrical scat- 
terers randomly placed inside a cylinder or radius R centered 
at position ro where the source is lying in the direct situation 
(surrounded by an exclusion surface of radius Ro). Detection 
wth a single point sensor lying at position in the far field 
[subfigure (a)] or a closed cavity (TRM) [subfigure (b)]. Real 
scatterers positions are displayed on that figure. 



electromagnetic problem is scalar. The source lies at the 
center of the cluster denoted by ro and is surrounded 
by an exclusion domain with radius Ro = 10 nm, small 
enough to preserve near-field interactions between the 
source and the surrounding scatterers. A single point 
antenna lies in the far field, outside of the system, at 
position = (4, 4) urn. To solve Maxwell's equations 
and to compute the electric field at the position of 
the atenna, we use the dyadic Green function G that 
connects the electric dipole at position ro to the radi- 
ated electric field at position Yd through the relation 
E(rd,u>) = noLu 2 G(rd,ro,Lj)po. To proceed, we perform 
a coupled-dipolc numerical computation. The field ex- 
citing scatterer number j is given by the contribution of 
the dipole source and of all other scatterers, leading to a 
set of TV self-consistent equations [20|, [2l[ : 



Ej = fi u) G (rj,r ,u))p + a(u})k 
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where Vj is the position of scatterer number j. 



Go (r,r» 



i/4 (k |r ■ 



r'|) is the free space Green 



function with H the Hankel function of first kind and 
order zero. This linear system is solved numerically for 
the configuration of Fig. Q] (a). Once the exciting electric 
field on each scatterer is known, it is possible to com- 
pute the scattered field at any position and in particular 
at the position of the antenna. In this numerical ap- 
proach, near-field and far field dipole-dipole interactions 
and multiple scattering are taken into account rigorously. 

In the time-reversed situation, the point source at po- 
sition ro is removed and the antenna becomes an active 
dipole source with an amplitude pd (cj) proportional to 
the complex conjugate of the recorded field. The same 
procedure as for the direct situation is used to compute 




FIG. 2: (Color online) Normalized broadband time-reversed 
intensity (A/T = 10) as a function of space for a single 
point antenna [subfigure (a)] and for a closed cavity [subfig- 
ure (b)]. Normalized monochromatic time-reversed intensity 
(A/r = 0) as a function of space for a single point antenna 
[subfigure (c)]. 

Figure [5] (a) displays the time-reversed intensity inte- 
grated over the frequency range [fi — A, il + A], defined 

as 



Itr (r) 



n+A 



Etr (r, uj) dw 



n-A 



(1) 



We clearly distinguish a focal spot centered at ro, the 
position of the source in the direct situation. The width 
of the focal peak can be estimated at A/6 (full width at 
half maximum) . This example demonstrates the feasibil- 
ity of subwavelength focusing in an open strongly scat- 
tering medium using a single point antenna. In order to 
determine the key parameters that explain the focusing 
process, we have performed numerical experiments. 

In a first step, we have replaced the single point an- 
tenna by an array of detectors placed in the far field 
forming a closed cavity around the system as shown in 
Fig. [T] (b), in analogy with the concept of TRM devel- 
oped for acoustics [7[. In such an ideal situation, the 
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time-reversed field is given by the imaginary part of the 
Green function of the scattering medium [22|: 



the disordered medium) 



E TR (r, u) = -2ifi uj 2 Q [G (r, r„, u)] p* 



(2) 



The corresponding time-reversed intensity, computed us- 
ing Eq. (JTJ is plotted in Fig. [5] (b). The pattern is sim- 
ilar to that in Fig. [2] (a) , with a width of the focal spot 
still on the order of A/6. We therefore conclude that 
time reversal at a single point antenna is here as efficient 
as time reversal in a closed cavity in terms of focusing 
performances. This is a consequence of strong multiple 
scattering in the disordered medium, that even in the 
case of an open system creates conditions for time rever- 
sal that are very close to that known for closed chaotic 
cavities 
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Nevertheless, strong multiple scattering and disorder 
are not enough to ensure focusing with super-resolution. 
In a second step, we have performed similar simulations 
with an exclusion domain of radius Rq = A around the 
focal point ro (not shown in this Letter for the sake of 
brevity) . This exclusion domain prevents near-field inter- 
actions between the source and the scatterers to occur. A 
time-reversed focal peak is still visible in these conditions, 
but with a diffraction limited width of A/ 2. This shows 
that near-field interactions with the surrounding scat- 
terers is necessary in order to generate a subwavelength 
concentration of light in the focus region. This mecha- 
nism has already been put forward in the context of time 
reversal of microwaves (l7j . It is interesting to note that 
the same near-field interactions are at the origin of the 
so-called Co correlation in speckle patter ns p roduce by a 
point source in a disordered medium 2J, |25| ) . Since this 
correlation is connected to the fluctuations of the local 



density of states [2a, [27[ , the latter being given by an ex- 
pression similar to Eq. ([2]) with r = ro, a close connection 
with time reversal focusing might be established (that is 
beyond the scope of this Letter). 

A crucial feature of time reversal is the spectral band- 
width. In a third step, we have studied the quality of 
focusing versus the bandwidth measured by the parame- 
ter A. In Fig. H](c), we show the time-reversed intensity 
/tr (r, fi) = |£tr (r, S!)| 2 for time reversal at a single 
point antenna and at the single frequency fl (which is 
equivalent to monochromatic phase conjugation). The 
signal-to-background ratio is very weak and no focal spot 
emerges in the intensity map. This demonstrates the cru- 
cial role of wideband time reversal for focusing in one re- 
alization of a disordered medium. A key question is the 
determination of the optimal bandwidth. Intuitively, it 
should be determined by the spectral correlation (or co- 
herence) bandwidth <5 co h of the speckle pattern produced 
in the disordered medium. Let us define the spectral 
field-field correlation at point ro as (the brackets (■ • ■) 
denote the average over an ensemble of configurations of 
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(3) 



and the coherence bandwidth S co h as the width of this 
correlation function. A self- averaging process can be ex- 
pected when A 3> (5 co h, giving a criterion for the band- 
width of the time reversal process ensuring focusing in 
one single realization of the medium. To compute C (5), 
we have generated numerically n — 112 random config- 
urations as that in Fig. [T] (a), and calculated the field 
at ro for plane-wave illumination. We have obtained a 
coherence bandwidth <5 co h ~ 0.1T, one hundred times 
smaller than the time-reversal bandwidth A = I0r. This 
amounts to manipulating A/<5 co h — 100 spectral degrees 
of freedom (compared to 1 in the monochromatic case) , 
explaining the improvement in the averaging process and 
the signal-to-background ratio. 

In order to illustrate the relevance of spectral correla- 
tions, we show in Fig. [3] maps of the time-reversed in- 
tensity for different frequencies around the reference fre- 
quency f2, and for a system in the multiple scattering 
regime (top row, optical thickness b ~ 16.8) or in the 
single scattering regime (bottom row, optical thickness 
b ~ 2). The three maps in the multiple scattering regime 
[subfigures (a), (b) and (c)] are clearly different, thus 
corresponding to three uncorrelated frequencies. Indeed, 
the coherence width in this case is <5 co h — 0.LT, small 
compared to the frequency steps between each map. In 
the single scattering regime [subfigures (d), (e) and (f)], 
the three maps present some similarities. The coherence 
width is <5 co h ~ 3r, larger than the difference in frequen- 
cies in the maps that remain correlated. This illustrates 
the positive role of multiple scattering in the efficiency of 
a broadband time-reversal process: a reduction of the co- 
herence width means a reduction of the bandwidth nec- 
essary to get self averaging in a single realization. An 
important conclusion of this work is that this behavior 
is preserved for focusing at subwavelength scale in the 
presence of near-field interactions and strong disorder. 

In summary, broadband time reversal at a single point 
antenna, in conjunction with near-field interactions and 
multiple scattering, is an efficient technique for spatial 
focusing of light at subwavelength scale in a disordered 
medium. The time reversal bandwidth necessary for self 
averaging has been connected to the spectral correlation 
width of the speckle pattern in a regime that was unex- 
plored so far. These results should give new perspectives 
for super-resolved optical imaging in complex media, and 
for the coherent control of single nanosources, including 
quantum emitters, or nanoscalc absorbers. 
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FIG. 3: (Color online) Normalized monochromatic time-reversed intensity Jtr (r, Q + 5). The first three subfigures correspond 
to the case of a system with a large optical thickness b = 2R/Ib (fi) = 16.8 [i.e. the system shown in figure [T] (a)] for three 
different frequencies given by 8/T = [subfigure (a)], 8/T = 0.5 [subfigure (b)] and 8/T — 1 [subfigure (c)]. The last three 
subfigures correspond to the case of a system with an optical thickness close to unity b = 2R/£b (fi) = 2.1 [i.e. a small section 
of the system shown in figure [1] (a)] for the same three different frequencies [respectively subfigures (d), (e) and (f)]. In all 
subfigures, the white dashed-dotted circle represents the boundary of the scattering cluster. The position of the single point 
antenna can be guessed in the top-right corner of the maps because of the presence of interferences fringes. 
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